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(57) Abstract: An improved method of acquiring seismic data using a plurality of vibratory seismic sources including the steps of 
deploying a seismic sensor; deploying a plurality of vibratory seismic sources at different source points; simultaneously actuating the 
seismic sources; acquiring seismic data attributable to the seismic sources using the seismic sensor; redeploying at least two seismic 
sources, one seismic source being thereby positioned at a source point previously occupied by the other seismic source; simultane- 
ously actuating the redeployed seismic sources; and acquiring seismic data attributable to the redeployed seismic sources using said 
seismic sensor. The present invention also involves an improved method of acquiring seismic data using a plurality of vibratory seis- 
mic sources, where each seismic source is capable of producing seismic energy within given frequency ranges, including the steps 
of deploying a seismic sensor; deploying a plurality of vibratory seismic sources at different source points; simultaneously actuating 
the seismic sources in such a manner that the frequency range of the seismic energy produced by one seismic source is substantially 
outside the frequency range of the seismic energy produced by another seismic source; and acquiring seismic data attributable to the 
seismic sources using the seismic sensor. 
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IMPROVED SEISMIC SURVEYING METHOD 

Background of the Invention 

This invention relates generally to seismic surveying methods, and more 
particularly to an improved seismic surveying method of using a plurality of 
vibratory seismic sources. 

Seismic vibrators have been used for many years on land to acquire 
seismic data and many companies have ongoing efforts to utilize similar sources 
in marine environments. The geophysical and environmental benefits of using 
these types of seismic sources are well known. 

When seismic data is acquired utilizing a plurality of vibratory seismic 
sources, the vibrators are conventionally organized as a travelling source array. 
The vibrators are typically placed around or along a source point (also referred to 
as a "vibrator point" or a "vib point") with a particular separation distance, such as 
40 meters. The vibrators then generate a certain number of sweeps that are 
received by a plurality of seismic sensors, recorded and stacked (i.e. combined) 
to produce a seismic data trace for each particular source point/receiver point 
pairing. The vibrators then travel as a group to the next source point where they 
are used in a similar manner. 

There are several known problems with acquiring seismic data using 
seismic vibrators, however, including the need to acquire large numbers of 
relatively-long records for each source point/receiver point pair to produce 
seismic data having a sufficiently high signal to noise ratio. Other known 
problems with seismic data acquisition using seismic vibrators include harmonics, 
ground coupling differences, baseplate flexures, and source array effects. 
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Efforts have been made to address these problems, and one promising 
approach has been the simultaneous use of multiple vibrators at different source 
points, with each vibrator producing separable, encoded sweeps. One method 
using this approach, referred to as the High Fidelity Vibroseis Source ("HFVS") 
method, has been developed by Mobil Oil Corporation and Atlantic Richfield 
Company and is described in U.S. Patent Nos. 5,550,786 (August 27, 1996); 
5,570,833 (December 30, 1997); 5,715,213 (February 3, 1998); and 5,721,710 
(February 24, 1998), all incorporated herein by reference. The HFVS method 
was developed primarily to improve the fidelity of vibroseis data. 

The HFVS method may be described, in principle, as comprising the 
following steps: 

1 . Measuring the vibrator motion S for each vibrator and each sweep, 
typically using an accelerometer mounted to the vibrator base-plate. 
The measured signal S is related to the true vibrator output U and a 
minimum phase transfer function T,. In the frequency domain, the 
equation describing the measured signal S is: S = U*Tl 

2. Recording the seismic data R. This seismic data represents the 
multiplication in the frequency domain between the earth reflectivity 
E, the vibrator output U and a minimum phase transfer function T 2 : 
R = U*T 2 *E. 

3. Obtaining the earth reflectivity at the vibrator location by multiplying 
the record R with the inverse of the vibrator motion U: R/U = Ti/T 2 *E 

For an array of 4 vibrators, V n , V 2j V 3 , and V 4l sweeping simultaneously, 
the geophone response R is described in the frequency domain by the following 
linear equation: R = mn*hi + mi 2 *h 2 + rrh 3 *h 3 + m 14 *h 4 . This equation contains 4 
unknowns, hi, h 2l h 3) and h 4 (the earth response at the vibrator positions V 1f V 2) 
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V 3( and V 4 ) and contains the known values R (the geophone response) and mn, 
nrh2, mn, and m 14 (the measured signals). 

The unknowns hi, h 2t h 3) and h 4 can be determined if another 3 sweeps are 
generated at the same locations and if the sweeps are encoded in such a way 
that the measured signal matrix is invertable. The system of linear equations is: 

Ri = mn*hi + mi2*h 2 + m 13 *h 3 + m 14 *h 4 
R2 = m 2 i*hi + m 22 *h 2 + m 23 *h 3 + m 24 *h 4 
R3 = m 31 *hi + m 32 *h 2 + m 33 *h 3 + m 34 *h 4 
R 4 = m 4 i*hi + m 42 *h 2 + m 43 *h 3 + nru^hu . 

In matrix notation, this can be written as: 

R = m x h 

where 

R- *" 

The typical implementation of the HFVS method in the field involves one 
array or group of vibrators, often four, spread out on an equal number of 
consecutive stations or source points. The vibrators sweep a certain number of 
sweeps, let say N (N being greater than or equal to the number of vibrators) at 
the same locations. The sweeps have the same frequency content but the phase 
is differently encoded to assure that the matrix M is invertible. After N sweeps, 
the vibrators move up a number of stations equal to the number of vibrators and 
repeat the sequence. 



m = 



m n m 12 m 13 m, 4 
m 2I m 22 m 24 
m 31 m 32 m 33 m M 
m 41 m 42 m 43 



and h = 



h 2 



3 



WO 00/72049 



PCT/IBOO/00654 



This implementation of the HFVS method has typically performed well in 
areas with shallow targets and good signal to noise ratios. For deeper targets or 
poor signal to noise areas, the standard implementation of the HFVS method may 
not perform well. The number of traces required for each source/receiver pair 
(the 'told") is also quite high, making acquisition of seismic data using this 
method relatively expensive. 

It is therefore desirable to implement an improved method of acquiring 
seismic data using a plurality of vibratory seismic sources that overcomes 
problems exhibited by prior art seismic data acquisition methods. 

An object of the present invention is to provide an improved method of 
acquiring seismic data using a plurality of vibratory seismic sources. 

An advantage of the present invention is that for the same acquisition 
effort and expense, seismic data having a higher signal to noise ratio may be 
obtained. 

Another advantage of the present method is that if coherent noise in the 
seismic data is band limited, it may be attenuated only in a particular frequency 
range, leaving the remaining frequency components of the seismic data 
unaffected. 

Summary of the Invention 

The present invention provides an improved method of seismic surveying 
using a plurality of vibratory seismic sources, the method including the steps of: 

deploying at least one seismic sensor; 

deploying a plurality of vibratory seismic sources at different source 

points; 



4 



WO 00/72049 



PCT/IB00/00654 



simultaneously actuating said seismic sources; 

acquiring seismic data attributable to said seismic sources using 
said seismic sensor; 

redeploying said seismic sources so that at least one of them is 
positioned at a source point previously occupied by another of them; 

simultaneously actuating said redeployed seismic sources; 

acquiring seismic data attributable to said redeployed seismic 
sources using said seismic sensor; 

decomposing said acquired seismic data into components 
attributable to each said seismic source; and 

stacking together components attributable to seismic sources 
located at a common source point. 

The invention and its benefits will be better understood with reference to 
the detailed description below and the accompanying drawings. 

Brief Description of the Drawings 

Figure 1 is a process flow chart showing steps associated with the 
inventive method; 

Figure 2 is schematic plan view of an exemplary seismic data equipment 
layout scheme; 
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Figure 3 is an exemplary amplitude versus frequency plot for a plurality of 
vibrators; 

Figure 4 is a fold distribution and seismic data section produced using a 
prior art seismic data acquisition method; 

Figure 5 is a seismic data section produced using the inventive station- 
shifting technique; and 

Figure 6 is a seismic data section produced using the inventive station- 
shifting and frequency splitting techniques. 

Detailed Description of the Invention 

Figure 1 is a process flowchart showing steps associated with the present 
method. The steps in Figure 1 will be discussed in connection with the schematic 
plan view of deployed land seismic acquisition equipment shown in Figure 2. 

As noted in Figure 1, typically the first step of the method sequence 10 is 
the "deploy seismic sensor" step 12. In Figure 2, a plurality of seismic sensors 
40, often geophones, are shown deployed along a plurality of parallel lines and 
connected to a data telemetry cable 42 which transmits the output of the sensors 
to a recording truck 44 where the acquired seismic data is recorded and often 
initially processed. Figure 2 depicts a typical 3D land seismic survey layout, but 
it is only one of a vast number of alternative seismic sensor deployment schemes 
that could be utilized in connection with the inventive method. 

Typically the second step of the method sequence 10 is the "deploy 
seismic sources" step 14. In Figure 2, source points 46 are represented as 
triangles and four vibrators V 2 , V 3l and V 4 (represented as circles) are being 
used in this seismic survey. Initially, in this example, vibrator Vi is located at 
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source point 48, vibrator V 2 is located at source point 50, vibrator V 3 is located at 
source point 52, and vibrator V 4 is located at source point 54. 

Typically the third step of the method sequence 10 is the "simultaneously 
actuate seismic sources" step 14, in which all four vibrators are simultaneously 
actuated to produce four successive sweeps each. In general, if there are N 
vibrators, each will be actuated to produce M successive sweeps, where M is not 
less than N. The vibrators will often be both phase and frequency encoded to 
provide enhanced signal separability. The phase encoding scheme for could, for 
instance, comprise the following: 



V! V 2 V 3 V 4 

Sweep 1 90 0 0 0 

Sweep 2 0 90 0 0 

Sweep 3 0 0 90 0 

Sweep 4 0 0 0 90 . 



Other methods for phase encoding the vibrator sweeps are described in 
the references incorporated earlier. This phase encoding assures the invertibility 
of the vibrator motion matrix. 

To enhance the separability of the signals, the sweep bandwidth is split 
among the vibrators using a frequency splitting technique to provide an additional 
degree of orthogonality to the source signals. If the sweep bandwidth required in 
a certain geologic area is between ^ and f 2 and an array of four vibrators are 
used to acquire the seismic data, the bandwidth may be split in the following way. 

Vi:[fi,fi + (f 2 -fi)/4]; 

V 2 :[fi + (f 2 -fi)/4 1 fi+(f 2 -fi)/2]; 

V 3 : [ fi + (f 2 - fi)/2, U + (f 2 - f 0*3/4 ]; and 

V 4 :[fi + (f 2 -fi)*3/4 ( f 2 ]. 
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Seismic vibrators are, however, typically hydraulically driven mechanical 
devices that lack the ability to rigidly cut off the production of seismic energy at 
any particular frequency. Generally they taper or ramp down the energy 
produced at the highest and lowest desired frequencies. To account for this 
behavior, referred to as the "sweep taper", a small overlap between the 
bandwidths for each vibrator may be desirable. This is shown in graphical form 
in Figure 3. 

As discussed above, the desired range of frequencies is divided by the 
number of vibrators (in this case four) and, in addition, each of the vibrators may 
be assigned different (and slightly overlapping) frequency ranges. In Figure 3, the 
desired range of frequencies is from 8 to 97 Hz. By following the frequency 
separation scheme described above, vibrator attempts to produce seismic 
energy matching the first curve 60 (8-31 Hz), vibrator V 2 attempts to produce 
seismic energy matching the second curve 62 (30-53 Hz), vibrator V 3 attempts to 
produce seismic energy matching the third curve 64 (52-75 Hz), and vibrator V 4 
attempts to produce seismic energy matching the fourth curve 66 (74-97 Hz). 
The first sweep taper region 68 overlaps the second sweep taper region 70 
between 30 and 31 Hz. Similar sweep taper overlaps occur between 52 and 53 
Hz and between 74 and 75 Hz. 

In terms of separability of the received signals, it would be preferable to 
eliminate any overlap of these taper zones. In many cases, however, it is desired 
to produce a source signal that is spectrally flat, i.e. that has substantially the 
same amplitude over the entire frequency range. If the signal is phase (as well 
as frequency) encoded, the spectral flatness benefits may overweight the 
reduction in separability this minor overlapping of the frequency spectra 
produces. The inventive method does not require that a spectrally flat signal be 
produced, however. In some cases, for instance, it may be preferable to weight 
or concentrate the seismic energy with respect to a particular frequency range or 
ranges, particularly if geologic conditions in a particular area substantially 
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attenuate reflected seismic energy outside a certain narrow frequency range or 
ranges. 

For each of the four sweeps produced when the seismic sources Vi to V 4 
are simultaneously actuated, the seismic sensors 40 will receive the seismic data 
in the "acquire seismic data" step 18. The seismic data is typically transmitted 
back to a recording truck 44 where it is recorded, preferably along with signals 
from respective accelerometers on each vibrator representative of the vibrator 
motion. 

When the required number of records have been obtained, the seismic 
sources are redeployed in the "shift seismic sources" step 20. Using one 
implementation of the station-shifting technique, vibrator V 4 is shifted to source 
point 56, vibrator V 3 is shifted to source point 54, vibrator V 2 is shifted to source 
point 52, and vibrator Vi is shifted to source point 50. The "simultaneously 
actuate seismic sources" step 22 and the "acquire seismic data" step 24 are then 
repeated the required number of times for that particular equipment layout. 

Typically, the vibrators will continue to produce seismic energy in their 
assigned split frequency ranges, as discussed above, but other frequency range 
assignment schemes can easily be envisioned. For source separability reasons, 
it is important that the frequency range of the seismic energy produced by one 
simultaneously actuated seismic source be substantially outside the frequency 
range of the seismic energy produced by another simultaneously actuated 
seismic source. Preferably at least half of the seismic energy produced by one 
seismic source occupies a different frequency spectrum than half of the seismic 
energy produced by another simultaneously actuated seismic source. As noted 
above, some overlapping of the frequency ranges may actually be desirable, but 
the purpose of this frequency splitting is to allow the received seismic data to be 
decomposed into components attributable to different source points based, at 
least in part, on their differing frequency contents. 
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Path 26 shows that this process is repeated until the required number of 
records have been acquired for each vibrator for each source point. The next 
time the seismic sources are shifted, vibrator V 4 is shifted to source point 58, 
vibrator V 3 is shifted to source point 56, vibrator V 2 is shifted to source point 54, 
and vibrator Vi is shifted to source point 52. After the required number of traces 
are acquired in this deployment setup, the vibrators are shifted one more station 
toward the top of Figure 2, and the required number of traces are acquired in this 
new deployment setup. After completing this fourth acquire seismic data step, it 
will be appreciated that each of the vibrators (V^ V 2 , V 3 , and V 4 ) have been 
deployed at source point 54. 

Each of the traces acquired by a particular seismic sensor while a vibrator 
was located at source point 54 will then be decomposed or separated in the 
"decompose acquired seismic data" step 28. Typically this decomposition 
process results in individual traces for each sweep each vibrator generated at the 
particular source point. If four vibrators are used and they each produced four 
sweeps at the source point, there would be 16 output traces, four each in four 
separate frequency band groups. Collectively, these four separate frequency 
band groups cover the entire seismic bandwidth required for the seismic survey. 

To increase the signal to noise ratio of the seismic data, these traces are 
combined ("vertically stacked") in the "stack decomposed seismic data" step 30. 
The traces may also be filtered prior to this stacking step. In some areas, noise, 
such as ground roll noise or power line noise, will be frequency band-limited. In 
these cases it may be desirable, for instance, to filter the frequency-limited traces 
containing this noise before they are stacked. This may be advantageous 
because noise attenuation methods may be used on the noise affected data 
without fearing that the noise attenuation methods may corrupt the seismic 
signals in the other sweep frequency ranges. This stacking process may consist 
of merely averaging all of the signals or more sophisticated stacking methods 
may be used such as the diversity stacking method discussed in U.S. Patent No. 
3,398,396 to P. Embree or the covariant stacking method discussed in our PCT 
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Patent Application PCT/GB98/0381 9 (WO 99/32903) to G. Baeten, both of which 
are incorporated herein by reference. It should be noted that the word "stacking" 
when used in this context simply refers to combining and does not imply normal 
moveout correction or reduction to zero offset. 

A significant advantage of the inventive method compared to prior art 
seismic data acquisition methods is that the signal to ambient noise ratio of the 
seismic data may be significantly improved. The signal to noise ratio of seismic 
data acquired using vibratory seismic sources may be calculated using the 
following equation: 

Signal/ambient noise = NV*SQRT(NS*L*W) 

where: 

NV = number of vibrators; 
NS = number of sweeps; 
L = sweep length; and 
W = sweep bandwidth. 

By increasing the number of vibrators deployed at each source point 
(albeit sequentially, not simultaneously), the number of sweeps and the sweep 
length can both be reduced while simultaneously improving the signal to noise 
ratio of the acquired seismic data. 

The separability of the sources is improved using the split bandwidth 
technique because the fundamental for each sweep is different and the first order 
harmonics generated by the third and fourth vibrator do not overlap the first and 
second sweep frequencies. The inventive method is capable of attenuating not 
only pneumatically introduced harmonics but also vibrator base-plate/earth 
interface introduced harmonics. 
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The amplitude spectrum produced by the vibrators also has larger values 
for narrow bandwidths and the same sweep length. 

Amplitude spectrum = Af*SQRT(T/4*W) 
where: 

At = amplitude of the fundamental; 
T = sweep length; and 
W = sweep bandwidth. 

Each of the vibrators is therefore able to transmit a greater quantity of 
seismic energy into the ground per time unit by limiting the bandwidth swept. 

Even if the sweep bandwidth used is the same for each vibrator (i.e. the 
frequency splitting technique is not used), the vertical stacking/station-shifting 
method will still improve the signal to ambient noise ratio of the acquired seismic 
data. When using the vertical stacking/station-shifting technique, any specific 
vibrator correlated noise will be attenuated because the seismic data associated 
with the particular source/receiver pair will consist of data associated with each of 
the vibrators. In addition, the data associated with any particular source/receiver 
pair will typically be acquired over a larger time window, thereby helping to 
attenuate any ambient noise that is time variant. 

The seismic data produced by the inventive method will then typically be 
subjected to other seismic data processing techniques such as filtering, 
migration, etc. that are well known in the seismic data processing art. 

The benefits associated with the inventive method were confirmed using 
both synthetic data and in a field experiment. The results of the field experiment, 
which was conducted in the Delaware Basin, Ward Country, West Texas, are 
shown in Figures 4, 5, and 6. It should be noted that in contrast to the 3D 



12 



WO 00/72049 



PCT/IB00/00654 



acquisition geometry shown in Figure 1 , these field experiments were conducted 
using a 2D acquisition geometry, where the sources (vibrators) and the receivers 
(geophones) are deployed along a common line. 

Figure 4 shows the results of a conventional 2D HFVS seismic survey. In 
this experiment, four vibrators were used, the sweep frequency was from 8 to 96 
Hertz, the sweep length was 10 seconds, the number of sweeps produced at 
each vibrator point was 8, the receiver interval spacing was 200 feet, the shot 
interval spacing was 50 feet, the listening time after each sweep was five 
seconds, and the number of stations each vibrator was shifted after completing its 
sweeps was 4 (200 feet). As can be seen in the fold bar chart 80 at the top of 
Figure 4, the fold produced by these acquisition parameters is approximately 60. 
The seismic data were then processed using a standard HFVS processing 
sequence (separation, spiking deconvolution, NMO correction, derive and apply 
statics, CMP sort, and stack) to produce the reference seismic data section 82 
shown in Figure 4. 

This reference seismic data section 82 can be compared against the 
improved seismic data section 86 shown in Figure 5. In this experiment, many of 
the parameters used were identical to those used in the previous experiment: four 
vibrators were used, the sweep frequency was from 8 to 96 Hertz, the receiver 
interval spacing was 200 feet, the shot interval spacing was 50 feet, and the 
listening time after each sweep was five seconds. The same vibrators, sensors, 
and recording equipment were used for this experiment and it was conducted 
over the same area as the first experiment. The sweep length was, however, 
reduced from 1 0 to 4 seconds, the number of sweeps at each vibrator point was 
reduced from 8 to 4, and the number of stations each vibrator was shifted after 
completing its sweeps was 1 (50 feet). As can be seen in the fold bar chart 84 at 
the top of Figure 5, the fold produced by these acquisition parameters was 
approximately 15 (a four-fold reduction from the previous experiment). The 
acquired seismic data were then processed using the same processing sequence 
to produce improved seismic data section 86. Many more coherent reflectors can 
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be seen in this improved seismic data section 86 than in the reference seismic 
data section 82. 

An experiment was also conducted using both the one-station station- 
- shifting technique and the frequency separation technique. The results of this 
experiment are presented in Figure 6. The only difference between this 
experiment and the previous experiment, the results of which are shown in Figure 
5, is that the four vibrators used separate sweep frequencies of 8-31, 30-53, 52- 
75, and 74-97 Hz, as discussed earlier. The fold was therefore unchanged, 
remaining at approximately 15 as shown in the fold bar chart 88 at the top of 
Figure 6. It should also be noted that the seismic data section 90 was created 
using onfy the 8-31 Hz bandwidth data because the frequency response in the 
test area was very narrow. The enhanced seismic data section 90 appears to 
provide an even clearer depiction of the subsurface geology than the improved 
seismic data section 86. 

It will be readily understood that the steps and processes associated with 
the disclosed embodiment of the present method are capable of a wide variety of 
alternative implementation methods and only a limited section from an actual 
seismic survey is discussed above. The described experiments imaged 
pressure-pressure transmission mode seismic energy, but the method is not 
limited to this particular seismic energy transmission mode and could, for 
instance, image pressure-shear converted transmission mode, shear-shear 
transmission mode, or multi-component seismic data. The seismic data acquired 
may be inverted using recorded vibrator output (using an HFVS-like method), 
may be inverted based on a theoretical or optimal vibrator output, or may be 
processed using other types of processing algorithms. The present method is 
also in no way limited or restricted to the particular order of steps described 
above. 
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CLAIMS 

1 . A method of seismic surveying using a plurality of vibratory seismic sources, 
the method including the steps of: 

deploying at least one seismic sensor; 

deploying a plurality of vibratory seismic sources at different source 

points; 

simultaneously actuating said seismic sources; 

acquiring seismic data attributable to said seismic sources using 
said seismic sensor; 

redeploying said seismic sources so that at least one of them is 
positioned at a source point previously occupied by another of them; 

simultaneously actuating said redeployed seismic sources; 

acquiring seismic data attributable to said redeployed seismic 
sources using said seismic sensor; 

decomposing said acquired seismic data into components 
attributable to each said seismic source; and 

stacking together components attributable to seismic sources 
located at a common source point. 

2. A method according to claim 1 , wherein said simultaneously actuating 
steps each comprise simultaneously actuating each seismic source M times at 
each source point, where M is not less than the number of seismic sources. 
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3. A method according to claim 1 or claim 2, further including the step of 
noise attenuating at least one of said components before said components are 
stacked together. 

4. A method according to any preceding claim, wherein the respective 
outputs of said seismic sources are recorded and used in processing said 
acquired seismic data. 

5. A method according to any preceding claim, wherein said seismic data is 
inverted using theoretical or optimal seismic source output. 

6. A method according to any preceding claim, wherein each said seismic 
source is capable of producing seismic energy within a respective frequency 
range and the frequency range of the seismic energy produced by one said 
seismic source is substantially outside the frequency range of seismic energy 
produced by another said seismic source when said seismic sources are 
simultaneously actuated. 

7. A method according to claim 6, wherein said seismic sources have sweep 
tapers and a sweep taper of one said seismic source overlaps a sweep taper of 
another said seismic source. 

8. A method according to claim 6 or claim 7, wherein the frequency range of 
one said seismic source has first order harmonics that do not overlap the 
frequency range of another said seismic source. 

9. A method according to any preceding claim, wherein said redeploying step 
comprising shifting said seismic sources one said source point in a common 
direction along a common path. 



16 



WO 00/72049 



PCT/IB00/00654 



1/5 



10 



14" 



18" 



Fig.1. 

Deploy Seismic 
Sensor 



Deploy Seismic 
Sources 



I 



Simultaneously 
Actuate Seismic 
Sources -* 



Acquire 
Seismic Data 



Shift Seismic 
Sources 



Simultaneously 
22-^ ^Actuate Seismic 
Sources 



Acquire 
Seismic Data 



28 



Decompose 

Acquired 
Seismic Data 



12 



-16 



~l 



20 



P 6 

I 



24 



I I 



Stack 
Decomposed 
Seismic Data 



30 



WO 00/72049 



PCT/IB00/00654 



2/5 



Fig.2. 



1 


40^ 40 


58 — -± 
56 J 


—-42 


V 4 AgL_-54 

nY — 52 


Kd — Q"^~44 


Vl — 48 




|~-46 





Fig.3. 




Freq. 



WO 00/72049 



3/5 



PCT/IB00/00654 



Fig.4. 



120 

ioo-~- 

2 80 

o 60 

40 



/ 



80 



120 



T>100_ 

"80 ^ 

Ul lll ll ll llll l HH ll Ul l Ull lll H l lllllllll I Mill II lllllilihiiiiiillll lllllllllliiiinim mil m — "^v^ .O 



20" 



Hi!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 



40 



U 1 1 1 M I II I 



-60 lu 
•20 



Col Loc 
0.0 



82 




Col Loc 
0.0 



WO 00/72049 



PCT/IBOO/00654 




WO 00/72049 



PCT/IB00/00654 




INTERNA 



10 



AL SEARCH REPORT 



al Application No 

PCT/IB 00/00654 



A. CLASSIRCATION OF SUBJECT MATTER 

IPC 7 G01V1/00 G01V1/02 



According to International Patent Classification (IPC) or to both national classification and IPC 



a FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 601V 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 

EPO-Internal , WPI Data, PAJ, INSPEC, COMPENDEX 



q DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ° Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



US 4 914 636 A (GARR0TTA ROBERT J) 

3 April 1990 (1990-04-03) 

abstract; claims 1,2,8,19; figures 2-4,6,7 

column 2, line 26 -column 4, line 34 

US 4 823 326 A (WARD ROGER M) 
18 April 1989 (1989-04-18) 
abstract; figures 2-5 
column 4, line 8 - line 17 
column 6, line 59 - line 66 
column 9, line 17 - line 29 



1,4,6,7, 
9 

2,3,5,8 



2 
8 



US 5 850 622 A (VASSILIOU ANTHONY A 
AL) 15 December 1998 (1998-12-15) 
abstract; figure 4 



-/- 



ET 



Further documents are listed in the continuation of box C. 



Patent family members are listed in annex. 



° Special categories of cited documents : 

'A' document defining the general state of the art which is not 
considered to be of particular relevance 

"E" earlier document but published on or after the international 
filing date 

"L" document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document pubfished prior to the international filing date but 
later than the priority date claimed 



T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the daimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art 

document member of the same patent family 



4 August 2000 


Date of mailing of the international search report 

11/08/2000 


Name and mailing address of the ISA 

European Patent Office, P.B. 561 8 Patentiaan 2 
NL - 2280 HV Rijswijk 
Tel. (+31-70) 340-5040, Tx. 31 651 eporrf. 
Fax: (431-70) 340-3016 


Authorized officer 

De Bekker, R 



page 1 of 2 



INTERNA 



AL SEARCH REPORT 



I Application No 

PCT/IB 00/00654 



C(Contlnuat)on) DOCUMENTS CONSIDERED TO BE RELEVANT 


Category ° 


Citation of document with indication, where appropriate, of the relevant passages 


Relevant to claim No. 


Y 


US 5 703 833 A (ALLEN KENNETH PAUL) 






30 December 1997 (1997-12-30) 






cited in the application 




A 
r\ 


OUOtl Otl) 1 1 yui CO x , £. 




Y 


US 3 885 225 A (ANSTEY NIGEL A ET AL) 


8 




20 May 1975 (1975-05-20) 




A 


i iyui c <j 


1,2,6 




column 6 line 11 — line 19 

\* V/ 1 Ullll 1 w j III IC XX III X J 




x 


W0 97 18488 A (MOBIL OIL CORP) 


1 

X 




22 Mav 1997 (1997-05-22) 






abstract* claims 1 3 

at/ o wi aw v* i a i mo x , ^ 




x 

A 


FR 2 338 499 A (BARBIER MAURICE) 

i i\ b j ju ~ j j i» v uni\u i mu i\ x w ^ / 






12 August 1977 (1977-08-12) 


A 


page 1, line 21 - line 38; figures 1,2 


2 




nanp 2 linp Q — linp ?fl 






nanp ? linp 35 — nan a 3 linp 1 

^sa^c ^ , ii i iv «j po^c «j j ii i icr x 






page 3, line 33 - line 39; claims 4,5 






page 5, line 15 - line 22 






(Jaytr v, I l lie &. x I I fie C/ 






narrp 1 linp 34 — linp *^ft 
payc x j ii lie? oh 1 1 lie jo 




A 


US 4 982 374 A (EDINGTON B LEON ET AL) 


1,2 




1 January 1991 (1991-01-01) 




paI limn 1 linp 1A — linp £1 
i»ui until ±, i i nc it i iiie u± 




A 


US 4 037 190 A (MARTIN I TNCOI N A^ 


Q 
O 




19 July 1977 (1977-07-19) 






abstract 






column 5, line 9 - line 22 




A 


US 5 721 710 A (SALLAS JOHN J ET AL) 


1 




24 February 1998 (1998-02-24) 






cited in the application 






abstract 











Form FCT/ISA/210 (continuation of second sheet) (July 1992) 



page 2 of 2 



INTERN. 



AL SEARCH REPORT 

n potent family members 



I Application No 

PCT/IB 00/00654 



Patent document 
cited in search report 


Publication 
date 


Patent family 
member(s) 


Publication 
date 




03-04-1990 


FR 


2622022 A 


21-04-1989 






AT 


65611 T 


15-08-1991 






AU 


2395388 A 


20-04-1989 






BR 


8805427 A 


20-06-1989 






CA 


1325840 A 


04-01-1994 






DE 


3863895 D 


29-08-1991 






EP 


0313459 A 


26-04-1989 






NX 


171345 B 


20-10-1993 



US 4823326 



18-04-1989 



NONE 



US 


5850622 


A 


15-12- 


-1998 


CA 


2242757 A 


14-05-1998 












EP 


0873527 A 


28-10-1998 












NO 


983130 A 


07-09-1998 












WO 


9820367 A 


14-05-1998 


us 


5703833 


A 


30-12- 


-1997 


CA 


2235411 A 


22-05-1997 












CN 


1201530 A 


09-12-1998 












EP 


0861450 A 


02-09-1998 












NO 


982146 A 


13-05-1998 












WO 


9718491 A 


22-05-1997 


us 


3885225 


A 


20-05- 


-1975 


GB 


1423366 A 


04-02-1976 












AU 


5825473 A 


23-01-1975 












CA 


984048 A 


17-02-1976 


wo 


9718488 


A 


22-05- 


-1997 


US 


5822269 A 


13-10-1998 












CA 


2233856 A 


22-05-1997 












rw 


L£\J£dDC A 


16-12-1998 












EP 


0861448 A 


02-09-1998 












NO 


982144 A 


13-05-1998 


FR 


2338499 


A 


12-08- 


■1977 


NONE 






US 


4982374 


A 


01-01- 


■1991 


DE 


69008998 D 


23-06-1994 












DE 


69008998 T 


01-09-1994 












EP 


0425250 A 


02-05-1991 


US 


4037190 


A 


19-07- 


1977 


US 


4058791 A 


15-11-1977 












US 


4189704 A 


19-02-1980 


US 


5721710 


A 


24-02- 


1998 


US 


5719821 A 


17-02-1998 












BR 


9706780 A 


04-01-2000 












CA 


2237408 A 


26-03-1998 












EP 


0927369 A 


07-07-1999 












NO 


982271 A 


17-07-1998 












WO 


9812578 A 


26-03-1998 



Form PCT/1SA/210 (patent tamiy annex) (July 1992) 



